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We present a model of polarization swings in blazars from axially symmetric blobs propa- 
gating on curved trajectories. If the minimum inclination of the velocity vector to the line 
of sight is smaller than T -1 , the polarization angle maximum rotation rate is simulta- 
neous with the polarization degree minimum and a spike in the total flux. By measuring 
the maximum rotation rate and the moment of the polarization maximum, it is possible 
to estimate the distance covered by the blob and thus its approximate position. We apply 
this model to the recent polarization event in blazar 3C 279. 
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1. Introduction 

Radio to optical emission of relativistic jets in Active Galactic Nuclei (AGNs) is lin- 
early polarized and this fact gives the strongest support for its synchrotron origin. 
Both polarization degree (PD) and polarization angle (PA) are strongly variable. 
Most of the time their behaviour is rather chaotic, but occasionally they show co- 
herent events. Two kinds of PA variations have been observed: systematic rotations 
of amplitude much larger than 180° on timescales of months to years (e. g. in PKS 
0721-115 1 and PKS 0521-365 2 ) and fast swings of amplitude close to 180° (e. g. in 
S5 0917+624 3 and BL Lac 4 ). 

Theoretical efforts to explain the PA rotation in AGN jets began with a model 
of accelerating clumps of matter 5 . A net transverse magnetic field component is re- 
quired to break the axial symmetry of the emitting region, so that the resulting PA 
rotation is smooth and gradual. It was later shown that other kinematic effects caus- 
ing variations of the viewing angle in comoving frame (jet bending or reorientation) 
could trigger a PA swing 6 . These models could explain PA swings of amplitude up 
to 180°. For larger PA rotations, models involving helical 7 or stochastic 8 magnetic 
fields have been proposed. 

In all these early works some kind of non-axisymmetric internal jet structure was 
assumed. Here we show, that it does not need to be the case. A symmetric emit- 
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ting region propagating on a bent trajectory, not necessarily helical, can produce a 
gradual PA swing. This is arguably the simplest way to reproduce several features 
in the simultaneous behaviour of PA and PD in some observational datasets with 
a minimum number of parameters. 



2. Planar bent trajectory 

Let us assume that for the time the PA is observed, the jet emission is dominated 
by a single compact emitting region ('blob') filled with magnetic fields axially sym- 
metric with respect to the velocity direction. It is known, that in such a case the 
polarization electric vector would be cither perpendicular or parallel to the pro- 
jected symmetry axis 9 , so it would be determined by the blob trajectory. To obtain 
a smooth PA rotation, the trajectory has to be curved and confined (at least locally) 
to a plane tilted to the line of sight by angle 9 m i n (see Fig. 1). The simplest possible 
case is that of planar curve of constant curvature radius R (circle) and constant 
blob velocity /3biob- When projected on the plane of the sky, it would appear as a 
curve bent by an angle equal to the observed PA rotation amplitude, i. e. ~ 180°. 
Interestingly, such large jet bends have been observed with VLBI on kpc scales in 
PKS 1510-089 10 and PKS 2136+141 11 . 

The PA rotation rate would achieve its maximum (dx/dt) b s ^ max at moment t 0} 
when velocity vector inclination to the line of sight 9 is minimal and equal to 8 m i n . 
For perpendicular magnetic fields, PD is highest for 6' = it/2 (as measured in blob 
co-moving frame), or sin 9 — Vj~ et (in external frame), where Fj ct is the Lorentz 
factor of the jet flow. PD would drop to for 9' = = 9, thus for # m ; n < rr t we 
expect a minimum of PD exactly at to- The behaviour of PD should be a symmetric 
function of time with two maxima at to ± Atriobs.max- 

By measuring (dx/dt) b s , max and Airiobs.max from the observational data, we 
can put the constraints on the jet parameters. Taking into account the light-travel 
effects, we obtain: 
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where f2 = (3biobc/R is the angular velocity of the blob velocity direction. By as- 
suming Tjet and Tbiob, we can calculate 9 m i n and R and find the distance covered 
by the blob between the PD minimum and maximum: 

Ar b iob = R x arccos ( J , et ) . (3) 
V cos 6» min / 



This gives us the lower bound for the position of the emitting region rbiob > Arbiob- 



August 28, 2012 0:21 WSPC/INSTRUCTION FILE Nalewajko2 



Polarization Swings from Curved Trajectories 3 




Fig. 1. The simplest case of blob trajectory that would produce a smooth PA swing of amplitude 
180°. Schematically shown is its orientation with respect to the line of sight (left) and the shape 
projected on the sky (right). The trajectory is confined to a plane inclined to the line of sight 
at angle m ; n . The 'Doppler cone' of opening angle ~ l/rj et is indicated. Black points on the 
trajectory mark the positions, at which the velocity vector is aligned with the Doppler cone, i. e. 
its inclination to the observer is 1/Fj ct in the external frame and 90° in the co-moving frame. At 
these points PD reaches its maxima. 

3. Application to a polarization swing in 3C 279 

Recently presented 12 polarimetric observations of blazar 3C 279 with KANATA and 
KVA telescopes show PA swing with reported amplitude of ~ 208° coinciding with 
a significant PD decrease to 3% (see Fig. 2). The PD minimum coincides well with 
the midpoint of PA rotation and a peak in the total flux at to = MJD 54888. Both 
before and after the swing PA shows little variability for a couple of months and 
is consistent with VLBI jet orientation Xjet ~ 60° 13 . to also marks a significant 
decrease in the optical and 7-ray flux. After t PD rises to achieve a maximum 
of 25% at ti = MJD 54908, during which both PA and the total optical flux are 
very stable. This behaviour was not parallelled before the swing, as PD was higher 
(~ 30%) and no clear peak can be seen. 
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Fig. 2. Left: V-band flux and polarimetric data for 3C 279 as observed with KANATA tele- 
scope. A model for polarization swing from the bent jet for observed maximum PA rotation 
rate ((ix/ett)obs.max = 20°/d, time elapsed between PD minimum and maximum Atn bs,max = 
t\ — to = 20 d and bulk Lorentz factor Pj ct = r^iob = 15 is shown with a solid line. Vertical 
dashed lines mark the moment of PD minimum to = MJD 54888 and the moment of PD max- 
imum t\ = MJD 54908. Horizontal dotted lines mark asymptotic PA values of —117° and 63°. 
Right: distance travelled by the blob between tg and t\, as calculated from Eq. (3) for rj et = r blob . 



We propose a scenario, in which there is a major bending in the jet at some rb on d- 
The high flux state was dominated by emission from the inner jet (r < rbend) and the 
flux state transition was caused by a sudden decrease in the central engine activity. 
Then the last portion of strongly emitting plasma travelled down the bent outer 
jet, producing a coherent PA swing with the successive PD maximum. We can thus 
use the bent-trajectory model to estimate rb en d- We adopt Ai no bs,max = <i — to = 
20 d and (rfx/^)obs,max = 20°/d. Assuming the Lorentz factor Tbiob — Tjct = 
15, consistent with VLBI measurements of the superluminal motions 14 , one can 
uniquely determine the orientation (# m j n — 0.7°) and curvature radius (R ~ 85 pc) 
of the blob trajectory and calculate the Stokes parameters and the Doppler factor 
as a function of observed time. We plot calculated polarization parameters and flux 
in Fig. 2. PA has been shifted by a constant value to obtain asymptotic values of 
— 117° and 63°. Note that the observational data matches the simulated curve with 
exception of small departures just before and after the swing. The observed PD 
has been scaled down by a constant factor to match the maximum value to the 
observed peak of ~ 28%. In this scaling the minimal value of <~ 2% at to roughly 
agrees with the observations. The curve shown on the flux plot is the Doppler 
boosting factor T> 3+a , where a = 1.75 is the optical spectral index {F v oc v~ a ). 
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It has been normalized by a constant factor to match the observed flux peak. We 
note, that if V ~ Tjet at t\ and T> < (1 + /3j c t)r.jct at to, the ratio of the Doppler 
factors is close to 2 and the ratio of the fluxes at to and t\ is about 25. Thus, the 
flux component from a single blob takes a form of a short spike. It matches well 
the observed feature at to, but cannot explain the long-term emission trend. We 
conclude that a rather small fraction of the observed flux has been emitted by a 
single blob. 

Using Eq. (3), we calculate the distance Arbiob covered by the blob between the 
PD minimum at t and the PD maximum at ti (see Fig. 2). Since 9 mln <C rr t , it 
scales approximately like Ar b iob ~ c Atnobs,maxPbiob • For r jot = Tbiob = 15, we 
expect the jet bend to be located at r bon d > Ar Mob ~ 5 pc. 

4. Discussion 

Our model is most suitable for PA swings of ~ 180° with simultaneous PD minimum 
and brightness spike, and with characteristic fast rise and slow decay of PD. We 
note that reported amplitude 208° of the polarization rotation in 3C 279 is based on 
the values measured immediately before the event, when there was a small rotation 
in the other direction. We also predict an apparent ~ 180° twist in a pc-scale jet, 
that in general could be verified with high-resolution VLBI observations. However, 
in the case of 3C 279 this is not plausible, since there is little variability in the mm 
wavelenghts and the optical/7-ray flare is produced within the mm photosphere. 

Larger observed amplitudes of PA rotation, e. g. ~ 240° in BL Lac 4 , could be ex- 
plained by helical trajectories. If the emission source is located in the Poynting-flux- 
dominated region, they would be determined by conservation of angular momentum 
inherited from the central engine. However, they could also arise in the interactions 
of the jet with non-uniform and variable environment on several-pc scales. There 
are two observational facts supporting this scenario in the case of 3C 279. First, the 
kpc-scale jet of continuously changes its positional angle and the apparent velocity 
of the radio knots on timescale of years 13 . Second, a PA rotation in the other di- 
rection has been observed in earlier epoch 15 , this challenges the models predicting 
fixed helicity sign. 

If the emitting region is large enough with respect to the curvature radius, time 
delays between the arrivals of photons emitted from different positions across the 
emitting volume would affect the variation of polarization parameters. This would 
provide a practical method to estimate the radius of the emitting region. We could 
also be able to place a constraint on the magnetic field distribution across the 
blob extension, in particular distinguish between globally ordered (toroidal) and 
disordered (chaotic) magnetic fields. 

Long and intensive polarimetric monitoring of bright blazars is very important, 
since polarization swings are fast and rare events, that cannot be predicted in ad- 
vance. A good coverage is essential for verifying competing theoretical models. In 
the case of 3C 279 an important question is, whether the observed polarization 



August 28, 2012 0:21 WSPC/INSTRUCTION FILE Nalewajko2 



6 Krzysztof Nalewajko 

swing was related to the transition from high to low 7-ray flux, or was it purely 
coincidental. We will only know this, when we observe another polarization swing 
in similar circumstances. 
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